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We explore the properties of magnetic excitations on the trellis lattice which is relevant for the so-called 
telephone-number compounds A14CU24O41 and the system SrCu2C>3. The trellis lattice consists of two-leg 
ladders which are coupled in a strongly frustrated fashion. We use the effective model obtained for a single 
two-leg spin ladder to calculate the two-dimensional one-triplon dispersion and the corresponding one-triplon 
contribution to the dynamical structure factor. Special attention is laid on signatures of the frustrating inter- 
ladder magnetic exchange. A detailed suggestion is made for an experimental detection of this exchange in 
inelastic neutron scattering experiments. 

PACS numbers: 74.25.Ha, 75.40.Gb, 75.10.Jm, 75.50.Ee 



Low dimensional quantum antiferromagnets display a va- 
riety of fascinating properties. The reduced dimensionality 
can either be realized by strongly anisotropic magnetic ex- 
changes or by a strongly frustrating topology of the system. 
The so-called telephone-number compounds A14CU24O41 and 
the system SrCu2C>3 are often discussed as realizations of ide- 
alized quasi one-dimensional two-leg ladders 1 ' 2 ' 3 ! 4 ' 5 ! 6 ! 7 . These 
one-dimensional structures build a two-dimensional trellis lat- 
tice (see Fig. 1), i.e. the coupling between the ladders is 
strongly frustrated. But the detailed influence of such a cou- 
pling and its quantitative strength are currently unknown. 



Quantum chemistry calculations have revealed a sizable 
value for the magnetic exchange between the one-dimensional 
ladder structures for the system SrC^O^. It is therefore an 
important issue to clarify the impact of the two-dimensionality 
of these systems on their magnetic properties. Theoretically, 
this question has been treated by mean-field theory using 
the bond-operator approach which is reliable when the rung- 
coupling J±_ dominates 9 . Exact diagonalization studies were 
performed on small clusters which provide only a limited 
number of points in momentum spaced. In this work we ap- 
proach this issue by studying the effects of the inter-ladder 
coupling J mt on the elementary excitations of the single lad- 
der and its dynamics. We concentrate on the paramagnetic 
case which is the experimentally relevant one. We give a de- 
tailed description for a possible experimental detection and 
determination of such a coupling by inelastic neutron scatter- 
ing. 




int 



FIG. 1: Heisenberg exchange couplings (J±,J\\,Ji n t) and the four- 
spin exchange coupling J cyc (circles) considered in the Hamiltonian 
Q}. The rung coupling (horizontal solid lines) is J±, the one along 
the legs of the ladders (vertical solid lines) is denoted by J11 and the 
coupling between adjacent ladders (dashed lines) is denoted by 7; nt . 



The couplings considered are shown in Fig. [T] They refer 



to the Hamiltonian 

H = Ho+H cyc +Hi n( 
Ho = J±_^Si,j-J_Si,f,R 

hi 

+^|| J! S,-J;tS/+1 j+1;t 

i,;;TG{L.R} 

#cyc = ^cycj^ C^ij^H-l J+I^l) (S(j;rSj+1 J+1;r) 

+ (S/J;LS i j;R)(S,' + lJ + l ;L S i+ lJ + l ;R ) 

~ (S,:;;LS,'+ 1 , /+ 1 ; r) (S; j;rS,- + 1 j + 1 ; l) 



#int — ^int^SijjR S,-+i + S; J-1;L 
'J 



(la) 



(lb) 



(lc) 
(Id) 



where the pair i, j 6 % denotes the rungs, not the sites, counted 
along the unit vectors e\ and <?2 shown in Fig. [2] The sub- 
scripts L and R stand for the left and right spin, respectively, 
on the particular rung. 

In order to explain the influence of various changes of pa- 
rameters we give below the formulae for the two-dimensional 
dispersion (£>k „ on a trellis lattice and for the dynamic struc- 
ture factor Sic,q(o)). The approach starts from the effective 
model that we have derived previously for a sing le ladder.^. 
It is given in terms of triplons, the elementary excitations of 
the spin ladder. In the following, we give a summarized de- 
scription how we derived the effective model for the single 
spin ladder. Technical details are given in Refs. ll3lll4 

We apply a particle-conserving perturbative continuous 
unitary transformation (CUT) which uses the states on iso- 
lated rungs as reference. The CUT maps the Hamiltonian of 
the single ladder to an effective Hamiltonian which conserves 
the number of triplons. In order to calculate spectral prop- 
erties and to treat the inter-ladder coupling the relevant 
observables have to be transformed by the same CUT. Here 
we use the perturbative realization of CUT which gives series 
expansions for the effective Hamiltonian and the effective ob- 
servables in terms of the expansion parameters x := Ju/J± and 
-*cyc : = J eye/ 'J jl m me thermodynamic limit. The obtained se- 
ries are extrapolated giving reliable results in a broad range of 
the parameter space 

For the present purpose, we need to know , the single- 
triplon dispersion along the ladder and the spectral weight a 2 k 
of the single triplon when it is excited by S z or another spin 
component. This spectral weight is known as function of the 
wave vector k. The coupling between the ladders is treated on 
the mean-field level starting from the ladders, not the isolated 
rung dimers 9 , as unperturbed syste m 17 ' 18 . We stress that the 
hardcore property of the triplons is taken into account along 
the ladders but not perpendicular to them. 

The rungs of the trellis lattice form a Bravais lattice so that 
there is exactly one mode of given magnetic quantum num- 
ber S z per point in the first Brillouin zone. But the lattice has 
a skewed unit cell due to the fact that adjacent ladders are 
shifted with respect to each other by half the distance c be- 
tween two rungs in one ladder (see Fig. 2). Conventionally, 
one would express the momentum in units of the reciprocal 




a 



FIG. 2: Unit cell of a trellis lattice spanned by the unit vectors e\ 
and e2- The distance between two adjacent rungs on the same ladder 
is c. The distance between the legs of a ladder is a, the distance 
perpendicular to the ladders between the centers (blue dots) of two 
rungs on adjacent ladders is a^- 



basis vectors e\ and e%. But in this system it is more conve- 
nient to describe the momentum by its component along the 
ladders (k, measured in units of 2n/c) and by its component 
perpendicular to the ladder (q, measured in units of 2n/(2ci2) 
, cf. Fig. [2]). Note that as a consequence the dispersion co(k, q) 
does not have the periodicities one would naively expect. The 
periodicity in k and q separately is 2, not unity. But shifting 
k — > k + 1 and q — > q + 1 together reproduces the spectrum as 
it has to be.— 

We use the same approach as in Refs. [l7lfT8l to treat the 



coupling H mt between the ladders, 
are transformed according to 



The observables 5 



avy,z.R 



a 



s 



(2) 



where the dots stand for normal-ordered quadratic and higher 
terms in the bosonic operators. In the following we neglect 
these terms corresponding to multi-triplon contributions. Note 
that the symmetry 5"' ff = — Sflg holds on this level of approx- 
imation. The Fourier transform squared of a$ yields the one- 
triplon spectral weight a?. In this approximation, the total 
Hamiltonian of the trellis lattice is quadratic in terms of the 
operators t x ' y,z ^ and t x,y ' z . Neglecting the hardcore property 
allows us to diagonalize the Hamiltonian by a Bogoliubov 
transformation which is justified as long as the coupling 7i nt 
is small. 
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The following two-dimensional dispersion Cfy „ 1S obtained 



(Ok.q = (0, 







1 — ^^a 2 cos(7Z:(7)cos(7ri:) 



coy 



The two-dimensional dynamic structure factor reads 

S ktq (C0) = a 2 kq 8{(0-(0 Lq ) 



4,q 



= 2sin z 



2< 



Ttq- — I a~ k - 

2a 2 J dtk.q 



(3) 

(4a) 
(4b) 



The sine factor stems from the interference of the excitation 
processes from the left and from the right spin in each rung. 
Recall that a single triplon has an odd parity on the isolated 
ladder with respect to reflections about the centerlin e 1 ^ 12 . 

The data for an isolated ladder is shown in Fig. [3] for 
the one-triplon dispersion of a single ladder and in Fig. [4] 
for the one-triplon weight. The couplings and energies are 
given in units of the rung coupling J±, that is x := J«/J±, 
Xcyc '■= Jcyc/J-L- This constitutes the input data for the two- 
dimensional calculation. 
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FIG. 3: Dispersion (of of a single triplon for an isolated ladder. Up- 
per panel: One-triplon dispersion for x cyc = and x = 1 (black 
curve) and x = 1.5 (cyan/grey curve). Lower panel: One-triplon 
dispersion for x cyc = 0.25 and x = 1 (black curve) and x = 1.5 
(cyan/grey curve). 
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FIG. 4: Weight a\ of a single triplon for an isolated ladder. Upper 
panel: One-triplon spectral weight for x cyc = and x = 1 (black 
curve) and x = 1.5 (cyan/grey curve). Lower panel: One-triplon 
spectral weight for x cyc = 0.25 and x = 1 (black curve) and x = 1.5 
(cyan/grey curve). 




FIG. 5: The influence of the inter-ladder coupling x mt on the dis- 
persion CDfc q of a single triplon on the trellis lattice as shown in 
Figs. |2] for x = 1 and x cyc = 0. Upper panel: Data sets for q = 1 
and x; nt = {— 0.4; — 0.2;0.0;0.2;0.4}. Lower panel: Data sets for 
q = 2 and x mt = {-0.4; -0.2; 0.0; 0.2; 0.4}. 



Note that at x = 1 .5 the extrapolation of the weight at small 
momentum spuriously leads to slightly negative values. In this 
parameter regime the three-triplon continuum overlaps with 
the one-triplon dispersion leading to a possible decay and a 
finite life-time of the elementary excitations. Therefore we 
put the one-triplon spectral weight in this range to zeroi 2 ^ We 
refrain from displaying data for even larger values of x be- 
cause there is a growing body of evidence that x S (1, 1.5) is 
the relevant range for cuprate ladder a 4,5,6,7,19 ' 20 even though 
early quantum chemical calculations 21 indicated even larger 
values of i w 2. Moreover, our approach is less reliable be- 
yond X K, 1.5. 

In the following, we will first discuss the generic properties 
of the two-dimensional spectra. To this end, we concentrate 



on the case x = 1 and x cyc = 0. But the features obtained are 
similar for all couplings in the paramagnetic phase. The in- 
fluence of the inter-ladder coupling Xj nt on the dispersion of a 
single triplon is shown in Fig. [5] For q = 1.5 the effect of the 
trellis structure vanishes completely due to destructive inter- 
ference. The largest effects can be seen for q = 1 and q = 2. 
(We refer to q = 1 .5 and q = 2 instead of q — 0.5 and q = 
although this does not make any difference in the framework 
of our theoretical model. But it eases the contact to the exper- 
imental investigations, cf. e.g. Ref. 7, where one keeps away 
from zero momentum where no magnetic scattering occurs.) 
Clearly, deviations from q = 1 .5 lead to a lifting of the reflec- 
tion symmetry about k = ±0.5 and to a difference between 
q = 1 and q = 2. On the level of our description, a ferromag- 
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netic coupling — |x; nt | at q = 1 is equivalent to an antiferromag- 
netic coupling |jq nt | aXq = 2 and vice versa. The one-triplon 
gap is shifted for an antiferromagnetic coupling from k = 0.5 
to a lower momentum for q = 2 and from k = 0.5 to a higher 
momentum for q = 1 . 




-0.5 

k [rlu] 

FIG. 6: The influence of the inter-ladder coupling X[ nt on the spectral 
weight a| of a single triplon on the trellis lattice as shown in Figs. [2] 
for x = 1 and x cyc = 0. Upper panel: Data sets for q = 1 and xj n t = 
{-0.4;-0.2;0.0;0.2;0.4}. Lower panel: Data sets for q = 2 and 
Stat = {-0.4;-0.2;0.0;0.2;0.4}. 

The physical interpretation of the destructive interference 
at half-integer values of q or of k relies on the geometry of 
the trellis lattice, see Fig. [2] Consider, for instance, two hop- 
ping processes from the two rungs, where the arrows e\ and 
e2 end, to the rung, from where both arrows start. For half- 
integer value of q these processes have opposite sign so that 
they cancel each other. Hence no effect of 7; nt occurs. The 
same is true for two hopping processes from the two rungs in 
the ladder in the middle in Fig. [2] to one of the other rungs. 
For half-integer value of k these processes have opposite sign 
so that they cancel each other. Again, no effect of /j nt oc- 
curs. This explains why /; n[ has to be significantly large in 
order to imply a sizeable effect if at least one of the momenta 
q and k is close to a half-integer values. Additionally, Fig. 
[4] shows that for values of k which deviate substantially from 
half-integer values the spectral weight of the triplons is very 
low so that the inter-ladder coupling has only very limited in- 
fluence on their dynamics, even away from half-integer values 
of momentum k. 

Besides the energy of the elementary triplon excitation, also 
the spectral weight is affected by the inter-ladder coupling. 
The correspoding results are shown in Fig. [6] In accordance to 
the one-triplon gap, the maximum of the one-triplon spectral 
weight is shifted for an antiferromagnetic coupling from k = 
0.5 to a lower momentum k for q = 2 and from k = 0.5 to a 
higher momentum k for q = 1 . But the overall effect on this 
quantity is rather small due to the frustration. 

Next we discuss a possible experimental detection of the 
inter-ladder coupling in inelastic neutron scattering experi- 
ments. From the above, it becomes clear that measurements 
collecting events in the interval q G [1,2] do not allow to es- 



tablish the inter-ladder coupling. We consider the possible 
detection of an inter-ladder coupling to be an essential ques- 
tion because it concerns the dimensionality of the magnetic 
systems A14CU24O41 and SrCu2C»3. Therefore, we suggest to 
measure around q = 1 and around q = 2. The comparison of 
the two constant energy scans of k provide the data to deter- 
mine the value of Ji nt . We point out that the weight factor 

sin 2 (^qj^j is reduced only by 25% if one deviates from the 

optimum value of q = 1.5 (where the sine is unity) to q = 1 
or q = 2. In this argument, we assume the approximate ratio 
«2 = 1.5a. 
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FIG. 7: Constant ©-scan of the dynamical structure factor at ft) = 
0.257x w i m x = 1.5 and x cyc = 0.25 as a function of the momentum k 
in ladder direction. The black curves correspond to x m = 0.4 and the 
cyan/grey curves correspond to xj nt = 0.2. The uppper panel shows 
q = 1 and the lower panel shows q = 2. A Lorentzian broadening of 
7] = 0.05 is used in all panels. 



In our opinion, the best experimental way to determine the 
inter-ladder coupling by inelastic neutron scattering is to per- 
form constant energy scans S(k,q = ±1, (O = A+ 8) at mo- 
menta q = 1 and q = 2 for energies A + 8 slightly above the 
one-triplon gap A, i.e. 8 > 0. The resulting spectra will show 
characteristic asymmetries with respect to reflections about 
£ = 0.5, see Fig. 13 

The existing experimental data of magnetic excitations in 
A14CU24O41 and SrCu2C<3 are consistent with coupling con- 
stants x = 1.2-1.5 and x cyc = 0.2 - 0.25^&i In the fol- 
lowing, we discuss the parameter set x = 1.5 and x cyc = 0.25 
which has been deduced from recent inelastic neutron scatter- 
ing data for La4SrioCu24C»4i— . Typical curves are shown in 
Fig. [7] Note again that in our mean-field treatment a change 
from antiferromagnetic to ferromagnetic inter-ladder coupling 
\xint\ — > — \xmt\ leads to the same curves as in Fig.[7]where the 
momentum q is changed to q + 1 . 

The large anisotropy x = J\\/J± = 1.5 is questioned by 
quantum chemistry calculations^. It is argued that the mag- 
netic exchanges along the legs /|| and the rungs J± are al- 
most equal but a strong ferromagnetic inter-ladder coupling 
x int ~ —0.2 is present (the calculations were done for the com- 
pound SrCu2C<3). We want to point out that this scenario of an 
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isotropic exchange x w 1 is not consistent with our results. As 
discussed above, the effects of the inter-ladder coupling are 
washed out if one averages over the momentum q perpendicu- 
lar to the ladder direction. This has been done when fitting the 
experimental data in Ref. 7. Therefore, the deduced parame- 
ters x = 1.5 and x cyc = 0.25 will not change if one takes into 
account a finite coupling x- mt . 

Independently from the intra-ladder anisotropy x, the strong 
ferromagnetic inter-ladder coupling deduced in Ref.|||is an in- 
teresting new aspect of the physics in these materials. It will 
be interesting to see if the analysis of inelastic neutron scat- 
tering data using the results obtained in this work will yield 
values of the magnetic exchange J mt which are consistent with 
the quantum chemistry calculations. 

In summary, we have discussed the properties of magnetic 
excitations on the frustrated trellis lattice in the paramagnetic 
phase by calculating the one-triplon dispersion and the cor- 
responding one-triplon contribution to the dynamical struc- 
ture factor. The trellis lattice can be viewed as being built 
from coupling components of one-dimensional two-leg lad- 
ders. Technically, an effective model for a single two-leg lad- 
der is used as the starting point to discuss the influence of 
the inter-ladder coupling. The strong frustration of the lattice 
causes the system to be effectively one-dimensional, i.e. the 
magnetic excitation spectrum reveals many properties of a sin- 
gle two-leg spin ladder. The leading order of the inter-ladder 
coupling 7i nt interferes destructively around the minima of the 
dispersion in the ladder. Hence, the effects of small or even 
sizable inter-ladder couplings are generically small. Yet, they 
lead to typical asymmetries of the dispersion and the spectral 



weight of a single mode. 

Although the mobility of triplons is strongly reduced by 
frustration, we expect that two- or multi-triplon properties like 
the interaction between the triplons are affected in a stronger 
fashion. While the hopping of triplons cancels due to frustra- 
tion their interaction as induced by different couplings adds. 
Thus, it is an important challenge for future work to clarify 
the influence of the inter-ladder coupling on optical properties 
and on all other probes of multi-triplon properties. 

We have made a detailed prediction how to detect the mag- 
netic inter-ladder coupling in experimental systems like the 
so-called telephone-number compounds A14CU24O41 and the 
system SrCu?03. We have provided arguments that exper- 
imental averages over the momentum perpendicular to the 
ladder direction do not yield information about the inter- 
ladder coupling. Therefore, the recently deduced parameter 
set x = 1.5 and x cyc = 0.25 will not be affected by a finite 
inter-ladder coupling 7 . This experimentally established dif- 
ference between /|| and Jj_ calls for further quantum chemical 
analyses of the superexchange couplings. 

The study of the inelastic neutron scattering spectra at q = 1 
and q = 2 should clarify the sign and the size of the two- 
dimensional inter-ladder coupling in the experimental sys- 
tems. Thereby, the imporant issue of the dimensionality of 
these magnetic systems will be elucidated. 

We thank B. Buchner, A. Gossling, B. Lake, A. Lauchli, 
S. Notbohm, and A. Tennant for stimulating and helpful dis- 
cussions. 
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Recall that the overall total spectral weight at small momentum is 
small anyway. Additionally, the three-triplon contribution to the 
dynamical structure factor behaves like CO 3 close to the lower band 



edge leading to a long life-time and therefore a sharp resonance 
at small momentum. This high power results from the hardcore 
property of the triplons. 



